Cosmographic constraints on a class of Palatini f(R) gravity 
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Modified gravity, known as f(R) gravity, has presently been applied to Cosmology as a realistic 
alternative to dark energy. For this kind of gravity the expansion of the Universe may accelerate 
while containing only baryonic and cold dark matter. The aim of the present investigation is to 
place cosmographic constraints on the class of theories of the form f(R) = R — a/R n within the 
Palatini approach. Although extensively discussed in recent literature and confronted with several 
observational data sets, cosmological tests are indeed inconclusive about the true signal of n in this 
class of theories. This is particularly important to define which kind of corrections (infra-red or 
high-energy) to general relativity this class of theory indeed represent. We shed some light on this 
question by examining the evolution of the deceleration parameter q(z) for these theories. We find 
that for a large range of a, models based on f(R) = R — a/R" gravity in th e Palatini approach 
can only have positive values for n, placing thus a broad restriction on this class of gravity. 

PACS numbers: 



I. INTRODUCTION 

f(R) gravity examines the possibility of modifying Ein- 
stein's general relativity (GR) by adding terms propor- 
tional to powers of the Ricci scalar R to the Einstein- 
Hilbcrt Lagrangian [l[ . The cosmological interest in these 
theories comes from the fact that they can naturally 
exhibit an accelerating expansion phase of the universe 
without introducing dark energy. However, the freedom 
in the choice of different functional forms of f(R) gives 
rise to the problem of how to constrain the many possible 
f(R) gravity theories on theoretical and/or observational 
grounds. In this regard, much efforts have been devel- 
oped so far, mainly from the theoretical viewpoint 
(see also Refs. [|| for recent reviews). General princi- 
ples such as the so-called energy conditions [||, nonlocal 
causal structure have also been taken into account 
in order to clarify its subtleties. More recently, obser- 
vational cons traints from several cosmological data sets 
have been explored for testing the viability of these the- 
ories 

An important aspect that is worth emphasizing con- 
cerns the two different variational approaches that may 
be followed when one works with /(i?)-gravity theories, 
namely, the metric and the Palatini formalisms (see, e.g., 
Q). In the metric formalism the connections are as- 
sumed to be the Christoffel symbols and the variation 
of the action is taken with respect to the metric, whereas 
in the Palatini variational approach the metric and the 
affine connections are treated as independent fields and 
the variation is taken with respect to both. In fact, these 
approaches are equivalents only in the context of GR, 
i.e., in the case of linear Hilbcrt action; for a general 
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f(R) term in the action, they provide completely differ- 
ent theories, with almost no similarities and very distinct 
equations of motion 1 . 

In this brief report we will restrict ourselves to 
the Palatini formalism. By discussing the decelera- 
tion/acceleration history of the Universe through the 
evolution of the deceleration parameter q(z), we place 
stringent bounds on a class of f{R) theory given by 
f{R) = f{R) = R—a/R n . In particular, we find that the 
range of negative values of n predicts a behavior that is 
grossly inconsistent with cosmological observations (e.g. 
CMB [l5| ) , being therefore ruled out even if it is compat- 
ible with supernova, lensing statistics and local gravity 
constraints. 



II. PALATINI f(R) COSMOLOGIES 

f(R) cosmologies are based in the modified Einstein 
equations of motion derived from the dubbed f(R) grav- 
ity. The action that defines an f(R) gravity is given by 

S = 2^ / dix ^~9f{R) + Srn , (1) 

where k 2 = 8irG, g is the determinant of the metric ten- 
sor and S m is the standard action for the matter fields. 
Treating the metric and the connection as completely in- 
dependent fields, variation of this action with respect to 



These differences also extend to the observational aspects of these 
approaches. For instance, we note that cosmological models 
based on a power-law function all form in the metric approach 
fail in reproducing the standard matter-dominated era followed 
by an acceleration phase [13 . whereas in the Palatini approach, 
analysis of a dynamical autonomous systems for the same La- 
grangian density have shown that such theories admit the three 
post inflationary phases of the standard cosmological model @ . 
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the metric provides the field equations 

/ 



K*T U 



(2) 



while variation with respect to the connection gives 

Vp{f'V=9!r)=0. (3) 

In (J2| T„ v is the matter energy-momentum tensor which, 
for a perfect-fluid, is given by T^ v = (p m + p m )u^u u + 
Vm9p,v, where p m is the energy density, p m is the fluid 
pressure and w M is the fluid four-velocity. Here, we adopt 
the notation /' = df/dR and /" = d 2 f/dR 2 . The sec- 
ond equation §3§ give us the connections r^„, which are 
related with the Christoffel symbols of the metric 

SV by 



{%} + ^JJ {^du + Kd, 9^9 pa d a ) f . (4) 



Note that in ©, R^ u must be calculated in the usual 
way, i.e., in terms of the independent connection T^ u , 
given by (|4]) and its derivatives. 

We assume a homogeneous and isotropic Friedmann- 
Lemaitre-Robcrtson- Walker universe whose metric is 
9ta> = diag(—l,a 2 ,a 2 7 a 2 ), where a(t) is the cosmolog- 
ical scale factor. The generalized Friedmann equation, 
obtained from ([5]), can be written in terms of the red- 
shift parameter z = ao/a — 1 and the density parameter 

*lmn = KPrn.n I ( 3i?n ) aS 
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where p mo is the matter density today. In terms of these 
quantities, the trace of Eq. ^ provides an important 
relation: 



Rf -2f = -3H 2 a mo (l + zf 



III. COSMOGRAPHIC CONSTRAINTS 



(6) 



Cosmographic parameters linked to geometry play a 
key role in cosmology. In particular, the deceleration 
parameter, defined as q(a) = — ad/ a 2 , can be written in 
terms of H(z) and its derivative as 



1 dH 

H(z)~dz 



(l + z)-l 



(7) 



In what follows we use ([5J) and J7J) to place constraints on 
the free parameters (a, n) of one of the most extensively 
discussed f(R) theory, i.e., 



f(R) = R- — 
R n 



(8) 



Note that positive values of n provide IR modifications 
of GR while negative ones correspond to high-energy cor- 
rections. For n — this kind of gravity reduces to Ein- 
stein's general relativity plus a cosmological constant, 



Test 


Ref. 


n 


a 


SNe la {Gold) a 


f6] 


[-0.3,0.3] 


[-6.8,-2.2] 


SNe la (SNLS) a 


17} 


[-0.23,0.42] 


[2.73,10.6] 


SNe la (Union)" 


18} 


[-0.3,0.1] 


[1.3,5.5] 


H(zT 


M 


[-0.25,0.35] 


[2.3,7.1] 


Strong lensing 


[10J 


[-0.202,0.078] 


[2.89,4.67] 


SNLS + WMAP + SDSS 6 


tiJJ 


[-5.98,-2.12] xlO -6 




SDSS C 


[12J 


[-2.321,1.329] xlO" 5 





a + CMB - 
b + BAO 
c + CMB - 



BAO 

- matter power spectrum 
matter power spectrum 



TT CMB spectrum 



TABLE I: Current constraints on n and a for a f(R) gravity 
given by ((8]). The ACDM model corresponds to n = and 
a = 4.38). 



f(R) = R — 2A. It is worth noting that, for a theory 
like ([5]), Eq. © evaluated at z — imposes the following 
relation between fl mo and the pair (n, a) 



R 



n+l 



1 - 



(9) 



where Rq is the value of the Ricci scalar today. Note also 
that Rq can be determined from the algebraic equation 
resulting from equating ([5]) and ((6|) for z = 0. Hence, 
specifying the values of two of these parameters the third 
is automatically fixed. In other words, in the Palatini 
approach, one of the two free parameters n and a can be 
eliminated by the constraint (j9|). 

Recently, the observational viability of the functional 
form © have been discussed by several groups. Con- 
straints from distance measurements to type-la Super- 
nova (SNe la) [6j8[, H(z) determinations by differen- 
tial age method 3| , strong lensing statistics [13], mat- 
ter power spectra [ll| , Baryon Acoustic Oscillation peak 
(BAO) and CMB shift parameter have been derived 
providing the following intervals n £ [—0.3,0.42] and 
a e [1.3,10.6], at 99.7% (C.L.) (see Table 1). 

However, a glance through Table 1 shows that cosmo- 
logical tests are indeed inconclusive about the true sig- 
nal of n. This in turn gives rise to the following question: 
what kind of corrections to GR do theories like © indeed 
represent? The answer to this question surely depends of 
the signal of the parameter n. In what follows, we show 
that some light can be shed on this discussion just by 
examining the signal of q(z), as given by ([7]) and taking 
into account (JS])-©. 

In Fig. rjjwe show the evolution of the deceleration pa- 
rameter with z for several values of n and f2 mo = 0.26, 
which corresponds to the best-fit value found from a 
joint analysis involving SNe la + BAO + CMB + H(z) 
data [1, [3|. As can be seen from Fig. [U for positive val- 
ues of n, q(z) approaches the standard ACDM evolution 
at intermediary redshifts and the Universe has only one 
recent transition toward an accelerating universe. In this 
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FIG. 1: Cosmic deceleration/acceleration history as a func- 
tion of redshift for the f(R) gravity theory given by |(HJ. Sev- 
eral cases are shown for selected values of n and S7 mo = 0.26. 
While positive values of n provides a q(z) evolution very simi- 
lar to the standard ACDM model, negative values of n predict 
either an always accelerating universe or two periods of accel- 
eration separated by a deceleration phase. 



case, the transition redshift does not vary considerably 
with n. 

On the other hand, negative values of n have great in- 
fluence on the q(z) evolution. In particular, we find that 
f(R) models given by © with n < —0.255 predicts an al- 
ways accelerating universe, which is clearly incompatible 
with a past dark matter-dominated epoch (and the stan- 
dard i 2 / 3 law), whose existence is fundamental for the 
structure formation process to take place 2 . Note that for 
—0.255 < n < these modified scenarios predict two re- 
cent periods of acceleration separated by a deceleration 
phase. 

Fig. [5] shows the first transition redshift func- 
tion of n. We find that the first epoch of acceleration 
is pushed to higher redshifts as n approaches zero from 
negative values and that such an effect is extremely sen- 
sitive to this parameter. For instance, for n = —0.05 
the first transition happens at very high z (z t = 10 8 ) 
whereas for n < —0.08 it occurs during the matter epoch 
(z t < 10 3 ). It is also worth mentioning that, although we 
have plotted Figs. [T] and [2] for a single value of fi„ io (= 
0.26), we have indeed repeated our analysis for a larger 
interval (Sl mo e [0.1,0.3]). We found that the behavior 
of q(z) and z t is very similar to the one shown in both 
figures, which means that for the current allowed range 
of the matter density parameter, negative values of n 



2 This value of n is slightly dependent on the matter density pa- 
rameter in the [0.1, 0.3]. 



FIG. 2: The first transition redshift (z t ) as a function of n. 
Note that as n — > the epoch of the first acceleration phase 
is pushed to higher z. 



are incompatible with a past matter-dominated epoch. 
This means that to be consistent with cosm ological ob- 
servations from LSS and CMB data the only correction 
to GR from a /(i?)-gravity like (J5J) is a IR modification 
(n > 0). Note that an exception to this conclusion is the 
limit n — > (from negative values), which for all prac- 
tical purposes is equivalent to the late-time evolution of 
the standard ACDM scenario. 



IV. CONCLUDING REMARKS 

Cosmological models based on f(R) gravity may ex- 
hibit a natural acceleration mechanism without intro- 
ducing a dark energy component. However, in order to 
provide a realistic alternative describing the late-time dy- 
namics of the Universe, it is qualitatively well accepted 
(since the current observational data are inconclusive to 
this respect - see Table I) that such theories must provide 
only an IR correction to GR. In this paper, we have dis- 
cussed quantitatively this question by studying the evo- 
lution of the deceleration parameter with z and impos- 
ing the existence of a past matter-dominated epoch. For 
the class of f(R) gravity given by ©, we have shown 
that negative values of n (high-energy correction) pre- 
dict either an always accelerating universe (which is in- 
consistent with the past structure formation process) or 
two periods of acceleration separated by a deceleration 
phase. These results show that even passing the back- 
ground tests (such as supernova ) and escaping from local 
gravity constraints practically the entire branch of nega- 
tive values of n is ruled out by the cosmographic bounds 
discussed here. 



4 



Acknowledgments 

The authors acknowledge financial support from Con- 
sclho Nacional de Desenvolvimento Cienti'fico e Tec- 



nologico (CNPq - Brazil) and Instituto Nacional de Es- 
tudos do Espago (INEspago). NP and JS also thank 
FAPERN for the grants under which this work was car- 
ried out. 



[1] A. A. Starobinsky, Phys. Lett. B 91, 99 (1980); 
R. Kerner, Gen. Relativ. Gravit. 14, 453 (1982); 
J. D. Barrow and A. C. Ottewill, J. Phys. A: Math. Gen. 
16, 2757 (1983); J. D. Barrow and S. Cotsakis, Phys. 
Lett. B 214, 515 (1988); B. Li and J. D. Barrow, Phys. 
Rev. D 75, 084010 (2007). 

[2] M. Ferraris, M. Francaviglia and I. Volovich, Class. 
Quantum Grav. 11, 1505 (1994); A.D. Dolgov and M. 
Kawasaki, Phys. Lett. B 573, 1 (2003); E. E. Flana- 
gan, Phys. Rev. Lett. 92, 071101 (2004); G. J. Olmo, 
Phys. Rev. D 72, 083505 (2005); Phys. Rev. D 75, 
023511 (2007); Phys. Rev. Lett. 98, 061101 (2007); S. 
Capozziello, S. Nojiri, S.D. Odintsov and A. Troisi, Phys. 
Lett. B 639, 135 (2006); A. W. Brookfield, C. van de 
Bruck and L.M.H. Hall, Phys. Rev. D 74, 064028 (2006); 
T. Multamaki and I. Vilja, Phys. Rev. D 73, 024018 

(2006) ; N. J. Poplawski, Phys. Rev. D 74, 084032 (2006); 
A. de la Cruz-Dombriz and A. Dobado, Phys. Rev. D 74, 
087501 (2006); L. Amendola, R. Gannouji, D. Polarski 
and S. Tsujikawa, Phys. Rev. D 75, 083504 (2007); L. 
Amendola, D. Polarski and S. Tsujikawa, Int. J. Mod. 
Phys. D 16, 1555 (2007); W. Hu and I. Sawicki, Phys. 
Rev. D 76, 064004 (2007); Y.-S. Song, W. Hu and I. 
Sawicki, Phys. Rev. D 75, 044004 (2007); O. Berto- 
lami, C.G. Bohmer, T. Harko and F.S.N. Lobo, Phys. 
Rev. D 75, 104016 (2007); B. Li, J.D. Barrow and D.F. 
Mota, Phys. Rev. D 76, 104047 (2007); I. Navarro and 
K. Van Acoleyen, JCAP 0702, 022 (2007); S. Baghram, 
M. Farhang and S. Rahvar, Phys. Rev. D 75, 044024 

(2007) ; T. P. Sotiriou, Phys. Lett. B 645, 389 (2007); 
C.G. Bohmer, T. Harko and F.S.N. Lobo, JCAP 0803, 
024 (2008); S.A. Appleby and R.A. Battye, JCAP 0805, 
019 (2008); E. Barausse, T.P. Sotiriou and J.C. Miller, 
Class. Quantum Grav. 25, 105008 (2008); G. Cognola, 
E. Elizalde, S. Nojiri, S.D. Odintsov, L. Sebastiani and 
S.Zerbini, Phys. Rev. D 77, 046009 (2008); S. Nojiri and 
S.D. Odintsov, Phys. Rev. D 77, 026007 (2008); T. P. 
Sotiriou, S. Liberati and V. Faraoni, Int. J. Mod. Phys. 
D 17, 399 (2008); S. Nojiri and S.D. Odintsov, Phys. Rev. 
D 78, 046006 (2008); T.P. Sotiriou, Phys. Lett. B 664, 
225-228 (2008); C.S.J. Pun, Z. Kovacs and T. Harko, 
Phys. Rev. D 78, 024043 (2008); V. Faraoni, Phys. Lett. 
B 665, 135-138 (2008); S. DeDeo and D. Psaltis, Phys. 
Rev. D 78, 064013 (2 008); S.H. Pereira, C.H. Bessa and 
J.A.S. Lima, larXiv:0911.0622l [astro-ph.CO]; G. Cognola, 
E. Elizalde, S. Nojiri, S.D. Odintsov, P. Tretyakov and 
S.Zerbini, Phys. Rev. D 79, 044001 (2009); V. Faraoni, 
Phys. Rev. D 81, 044002 (2010); T. Harko, Phys. Rev. D 
81, 044021 (2010); T. Narikawa and K. Yamamoto, Phys. 
Rev. D 81, 043528 (2010); T. Multamaki, J. Vainio and I. 
Vilja, Phys. Rev. D 81, 064025 (2010); J. Kluson, Phys. 



Rev. D 81, 064028 (2010); E. Santos, Phys. Rev. D 81, 
064030 (2010). 

[3] S. Capozziello and M. Francaviglia, Gen. Relativ. Gravit. 
40, 357 (2008); T. P. Sotiriou and V. Faraoni, Rev. Mod. 
Phys. 82, 451 ( 2010); A. De Felice and S. Tsujikawa, 
larXiv: 1002.4 928 [gr-qc]. 

[4] J.H. Kung, Phys. Rev. D 53, 3017 (1996); S.E.P. Bergli- 
affa, Phys. Lett. B 642, 311 (2006); J. Santos, J.S. Al- 
caniz, M.J. Rebougas and F.C. Carvalho, Phys. Rev. D 
76, 083513 (2007); K. Atazadeh, A. Khaleghi, H. R. Sep- 
angi and Y. Tavakoli, Int. J. Mod. Phys. D 18, 1101 
(2009); O. Bertolami and M.C. Sequeira, Phys. Rev. D 
79, 104 010 (2009); J. Sa ntos, M.J. Reboucas and J.S. 
Alcaniz, i arXiv:0807.2443l [astro-ph] . 

[5] T. Clifton and J. D. Barrow, Phys. Rev. D 72, 123003 
(2005). M.J. Rebougas and J. Santos, Phys. Rev. D 80, 
063009 ( 2009); J. Santos, M.J. Rebougas and T.B.R.F. 
Oliveira, larXiv: 1004.25011 [astro-ph.CO]. 

[6] M. Amarzguioui, 0. ElgarOy, D.F. Mota and T. Mul- 
tamaki, Astron. Astrophys. 454, 707 (2006). 

[7] S. Fay, R. Tavakol and S. Tsujikawa, Phys. Rev. D 75, 
063509 (2007). 

[8] J. Santos, J. S. Alcaniz, F. C. Carvalho, N. Pires, Phys. 

Lett. B 669, 14 (2008). 
[9] F.C. Carvalho, E.M. Santos, J.S. Alcaniz and J. Santos, 

J. Cosmology Astroparticle Phys.: JCAP 09, 008 (2008); 
[10] X.-J. Yang and Da-M. Chen, Mon. Not. R. Astron. Soc. 

394, 1449 (2009). 
[11] B. Li, K. C. Chan and M.-C. Chu, Phys. Rev. D 76, 

024002 (2007). 
[12] T. Koivisto, Phys. Rev. D 73, 083517 (2006). 
[13] A. Borowiec, W. Godlowski and M. Szydlowski, Phys. 

Rev. D 74, 043502 (2006); B. Li and M.-C. Chu, Phys. 

Rev. D 74, 104010 (2006); M. Fairbairn and S. Rydbeck, 

JCAP 0712, 005 (2007); M. S. Movahed, S. Baghram and 

S. Rahvar, Phys. Rev. D 76, 044008 (2007); T. Koivisto, 

Phys. Rev. D 76, 043527 (2007); H. Oyaizu, M. Lima and 

W. Hu, Phys. Rev. D 78, 123524 (2008); C.-B. Li, Z.-Z. 

Liu and C.-G. Shao, Phys. Rev. D 79, 083536 (2009); 

K.W. Masui, F. Schmidt, Ue-L. Pen and P. McDonald, 

Phys. Rev. D 81, 062001 (2010). 
[14] L. Amendola, D. Polarski and S. Tsujikawa, Phys. Rev. 

Lett. 98, 131302 (2007). 
[15] E. Komatsu et al, Seven-year Wilkinson microwave 

anisotropy probe (WMAP) observations: Cosmological 

interpretation, arXiv:1001.4538 [astro-ph.CO]. 
[16] R. Jimenez and A. Loeb, Astrophys. J. 573, 37 (2002). 
[17] S.M. Carroll, V. Duvvuri, M. Trodden and M.S. Turner, 

Phys. Rev. D 70, 043528 (2004). 



